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Abstract. Let (M, g) be a two-dimensional compact boundaryless Riemannian man- 
ifold with Laplacian, A 9 . If e\ are the associated eigenfunctions of yj— A 9 so that 
—A g ex = A 2 e^, then it has been known for some time |18| that ||exllz, 4 (M) ~ A 1 / 8 , 
assuming that e\ is normalized to have L 2 -norm one. This result is sharp in the sense 
that it cannot be improved on the standard sphere because of highest weight spherical 
harmonics of degree k. On the other hand, we shall show that the average L 4 norm of 
the standard basis for the space "Hfc of spherical harmonics of degree k on S 2 merely 
grows like (logfc) 1 / 4 . We also sketch a proof that the average of Y^ 21 ^ 1 ||e.\||* 4 for 
a random orthonormal basis of "Hfe is O(l). We are not able to determine the max- 
imum of this quantity over all orthonormal bases of rlu or for orthonormal bases 
of eigenfunctions on other Riemannian manifolds. However, under the assumption 

that the periodic geodesies in (M, g) are of measure zero, we are able to show that 

1/8 

for any orthonormal basis of eigenfunctions we have that \\ e \j k lli 4 (M) = °(^j k ) ^ or 
a density one subsequence of eigenvalues Xj. . This assumption is generic and it is 
the one in the Duistermaat-Gullemin theorem [6] which gave related improvements 
for the error term in the sharp Weyl theorem. The proof of our result uses a recent 
estimate of the first author |20| that gives a necessary and sufficient condition that 
WexhHM) = oCA 1 / 8 ). 



1. Introduction. 

The purpose of this note is to introduce a new problem on LP norms of eigenfunctions 
on compact Riemannian manifolds (M,g). We prove some initial results on the problem, 
and also include some conjectures and heuristic remarks. 

The problem, roughly speaking, is to determine the asymptotic average of the L 4 
norms 1 1 e a. 1 1 4 of the elements of an orthonormal basis of eigenfunctions 

-A g e\ = \ 2 e\ 

of the associated Laplace-Beltrami operator. In practice it is simpler to consider the 
fourth power Weyl sums, 

(i.D ^ Y. Kiii 

where 

N(X) = #{A 3 < A}, 
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is the Weyl counting function. The asymptotics of depend on the entire orthonor- 
mal basis and, as will be seen below, can behave quite differently from the behavior of 
individual eigenfunctions in the basis. 

Before stating our results, let us recall the results on L p norms of individual eigen- 
functions. In 1988, one of us showed in [T8] that for 2 < q < oo and 

(1.2) a(q) = max( 2(1/2 - 1/q) - 1/2, ±(1/2 - 1/q) ) 

= f2(l/2- 1/2, q>6 

{ |(l/2 -1/q), 2<q<6, 

we have 

(1-3) IMU«(m) < A CT(9) , 

assuming as we shall do throughout that the eigenfunctions are L 2 -normalized so that 

ll e A|U 2 (A/) = 1, 

where the norms are taken with respect to the volume element, dV . This result is sharp 
since certain spherical harmonics on the sphere, S 2 , with the round metric saturate the 
estimate (| 1 .3|) . Specifically, when q > 6, L 2 -normalized zonal functions, Zk satisfy 

WZ.W^^k^-^-^, q>6, 

while the L 2 -normalized highest weight spherical harmonics, Qk — Ck{x\ + ix-i) k satisfy 

||Q fe |U^)«fc5 (1/2 - 1/9) , q>2. 

Both are eigenfunctions of the standard Laplacian on S 2 with eigenvalue A 2 = k(k + 1) 
in the above notation. Also, we are taking S 2 to be {(xi, x 2 , x 3 ) : x\ + x\ + x 2 = 1}, 
so that, as k —¥ oo, the Qk become highly concentrated on the equator where X3 = 0. 
The orthonormal basis of joint eigenfunctions of A g and of a^-axis rotations are generally 
denoted by Y/^, m = —k, . . . , k; in particular, Zk = Yq and Qk = Y k . 

Even though (| 1 . 2j) cannot be improved on the sphere, it is thought that for generic 
manifolds one has at least 

(1.4) \\ex\\mM) =°(A CT(9) ), as A — > 00, 

for a given q > 2. In [22] we showed that for generic (M,g) this is true for q > 6 (and 
also corresponding results for higher dimensions). This just followed from showing that 
under a certain generic condition on (M, g) one can improve the L°° estimate in (|1.3[) 
to be || e> 1 1 oo = o(A 1 / 2 ), which implies (|1.4j) for all q > 6 by interpolating with (|1.3j) for 
q = 6. The results in [2 2) were recently improved in |21j . The key point was to show that 
the bound ||ex||oo = (^(A 1 / 2 ) can only be obtained on (M, g) possessing a "peak point" or 
"pole" Zo with the property that a positive measure of directions in S* Q M exponentiate 
to geodesic loops which return to zq at some time. This behavior occurs at poles of a 
surface of revolution, since all meridians are closed geodesies through the pole, and in 
particular explains why the sup norm bounds are attained by zonal functions^ on the 
round sphere (see §3 below). 

Even though there are satisfactory results concerning (|1.4[) for relatively large expo- 
nents q > 6, much less is known for relatively small exponents 2 < q < 6. In this case, 
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it is thought that the enemy for (ll.4[) is maximal concentration along periodic geodesies, 
as occurs for the highest weight spherical harmonics. Using the formula for the Qk one 
checks that they have L 2 -mass bounded below on shrinking fc -1 / 2 neighborhoods of the 
equator 7 = {(xi, X2, x%) £ S 2 : x 3 — 0}. In [20] , (following an earlier result in [3]), the 
first author proved that for 2 < q < 6, (|1.4j) is valid if and only if this type of concen- 
tration does not occur. Specifically, a necessary and sufficient condition for (|1.4[) for this 
range of exponents is that 

(1.5) sup f \ex(y)\ 2 dV = o(l), 

■yen Jdist s (7,i/)<A- 1 /2 

where 77 is the space of all unit-length geodesies in M, and dist g ( •, • ) is the geodesic 
distance associated to the metric g. 

The goal of this paper is to show that even though on some manifolds there are 
eigenfunctions e\ having L 4 -norms of maximal size ~ A 1 / 8 as A — > 00, they are very 
sparse. Our first result of this type says that given any orthonormal basis {e^} of 
eigenfunctions with eigenvalues Ai < A2 < . . . on a two-dimensional compact Riemannian 
manifold (M, g) with a zero measure of periodic geodesies, one can find a density one 
subsequence of eigenvalues, {Aj fc }, for which 

(1-6) \\e x J LHM) =o(Xy k S ). 

By interpolation with the L 6 -estimate in (|1.3p and the trivial L 2 -estimate, this implies 
that we also have ||e^ = o(A°^) for every 2 < q < 6. Presently, we do not how 
to prove the corresponding results for q > 6, or how to obtain any results like this for 
higher dimensions n > 3. The condition that (M, g) have a zero set of periodic geodesies 
is generic and it is the assumption in the Duistermatt-Guillemin theorem [6], which 
involved a similar o-improvement of the error term in the Weyl formula. 

The assumptions that the periodic orbits are of measure zero of course is not valid for 
the sphere. Nonetheless, we can prove a much stronger result for the standard basis {Y^} 
on S 2 , even though, as we pointed out before, this eigen-basis has functions saturating 
(II. 3|) for each 2 < q < 00. 

To be more specific, we recall that the Laplace-Beltrami operator on S 2 C K 3 with 
the standard round metric has eigenvalues A 2 = fc(fc + 1) repeating with multiplicity 
2k + 1, meaning that the corresponding eigenspace T~Lk of spherical harmonics of degree 
k has this dimension. If we use longitudinal coordinates <p G [0,7r] and latitudinal ones 
9 € [0,27r] so that S 2 3 x = (sin cos 8, sin </>sin#, cos </>), then in these coordinates the 
standard basis for 'H.k has elements 

(1.7) Y*(cj>,6) - c mik P? (cos 0)e ime , -k < m < k, 

where PJJ 1 are Legendre functions and c m ^, are L 2 -normalizing constants. When m = 0, 
Yq is the zonal function and when m = ±fc it is a highest weight spherical harmonic 
of degree k. For this basis, we shall show that the average L 4 -norm is of size ~ (log fc) 1 / 4 , 
as fc — > 00, i.e., 

k 

(1-8) ^ £ / l^.logfc, fc>2, 
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which of course is much stronger than (jl.6p since it shows that there must be a density 
one sequence of eigenfunctions among this basis with L 4 -norms growing logarithmically 
with respect to the eigenvalues. It seems somewhat paradoxical at first that (|1.8p is valid 
for the standard basis on the sphere, while the same basis is the worst case for (|1.3|) . But 
this holds because the left side of (|1.8|) is a functional of an orthonormal basis rather than 
of individual eigenfunctions, and most elements have relatively small L 4 norms. It is 
doubtful that the {Y^} maximize this functional among orthonormal bases of spherical 
harmonics. In Section U we explain this further. 

These observations raise the following 

Problem Let dimM = 2. For which (M,g) (if any) does there exist an orthonormal 
basis of eigenfunctions for which there exists a positive density subsequence e\ jk so that 

\\ e *j k \\l 4 = ^(^jf 8 )-? Or is a result like (|1.6j) is valid on any compact surface? 

We prove (jTTSJ) by obtaining pointwise bounds for the £ 4 (m) norms of the basis elements 
of "Hfc. Specifically, we shall prove sharp estimates for 



By the inclusion £ 4 C £ 2 , this quantity is bounded by the corresponding £ 2 (m)-norm. On 
the sphere, the £ 2 (m) norm is independent of x, and, in fact, 



The £ 4 (m) norm is of this order of of magnitude for points x of distance 0(1/ k) from the 
poles where 4> = or ir, but in order to obtain (|1.8p much better estimates are needed. 
We shall obtain such an improvement, which turns out to be sharp, by using (jTTTJ) and 
well known asymptotics for the kernel of the projection onto the spherical harmonics of 
degree k,%k- Thus, we are very much using here special properties of S 2 . Our results 
can be thought of as a natural analog for S 2 of Zygmund's [27] theorem for the two-torus 
T 2 , which says that the eigenfunctions of its Laplace-Beltrami operator have uniformly 
bounded L 4 -norms. As we pointed out before, this is far from true on S 2 , but in an 
averaged sense it is almost true since the average L 4 -norms just grow like powers of logs 
of the eigenvalues. 

For general Riemannian manifolds of dimension n, the local Weyl formula says that if 
N is large enough and fixed then 



It would be interesting to see to what extent there is an improvement in the general case 
when one replaces this £ 2 -norm by £ q norms with q > 2 and to what extent results of this 
type perhaps depend on properties of the geodesic flow starting at a;. In a future work, 
we intend to carry out the analysis for round spheres of dimension n > 3 and certain 
surfaces of revolution. Understanding the case of general manifolds and to what extent 
these results might depend on x seems difficult at present. On the other hand, by using 




(1.9) 



\\Y*(x)\\ P{m) = V(2fc + l)/47r. 



(1.10) 
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estimates like (ll.3[) . one can see that for most points x € M, once can improve on the 
trivial consequence of (|1.10[) that 



1/9 



E M*)! 4 

|Aj-A|<JV 



< ^(»-l)/2_ 



For instance, if q = 4 and n = 2, then using (|1.3I) and Tchebyschev's inequality one sees 
that if C < oo is fixed then 



xeM:( E I^a, (a;)| 4 ) 1/4 > CA 1 / 2 
|A 3 -A|<JV 



which, not surprisingly, is exactly the size of the sets on which the highest weight spherical 
harmonics are concentrated. 

This paper is organized as follows. In the next section we shall present the proof of 
(|1.6j) . Then we shall turn our attention to the sphere S 2 and prove the much stronger 
bounds (JUSJ for S 2 . 



2. L 4 norms of generic eigenfunctions. 

In this section we shall establish (|1.6|) . Specifically, we shall prove the following 

Theorem 2.1. Let (M,g) be a two-dimensional compact Riemannian manifold. If 
S*M — > S* M is geodesic flow on the cosphere bundle, assume that the set 



(2.1) 



V = {(£,£) e S*M : $t(x,Q = (x,0, s ome t > 0} 



has measure zero in S*M with respect to the volume element. Then if e\. is an orthonor- 
mal basis of eigenfunctions, — Ae^ = Xj, with Ai < Aa < . . . there is a subsequence of 
eigenvalues Xj k satisfying 



(2.2) 

so that 
(2.3) 



lim 



#{A Jfc <A} 
N(X) 



1, 



To prove this we shall use an estimate from [3D] and arguments from [S] and [33] . The 
estimate from [3D] says given (M, g) as above there is a uniform constant C so that if 
-Ae A = X 2 e\ and N = 1,2,3, ... then 



\ex{x)\UV <CN-V 2 A 1 / 2 || eA ||J 



A I 



L 2 {M) 



sup 

ienJr x - 1/2 h) 



\ex(x)\ 2 dV 



L 2 (M)- 



Here dV = dV g is the volume element, 77 is the space of all unit-length geodesies, and 



Teh) = {2/SM : dist fl (y, 7 ) < s}, 
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denotes an e-tube about 7. By optimizing the choice of N, we see that the preceding 
inequality implies that 

(2.4) \\ex\\L H M) < CX^exfJ* sup ||e A ||$ r + C\\e x \\ LHM) . 

y ' yen y A 1/2 ' 

In addition to this we require the following result which is a simple consequence of the 
local Weyl law (see [lljh 

Lemma 2.2. Let (M,g) be a compact Riemannian manifold and let A 6 ^^(M) be a 
classical pseudo- differential operator on M of order zero. Then if Aq is the principal 
symbol of A, 



A,<A 



(2.5) V / \Ae Xj (x)\ 2 dV = (2ir)- n \ n / \A (x,£)\dVd£ + 0(A™ _1 ). 

' M JT*B 



Here, T*B C T*M is the ball bundle, {(a:,£) : J2 jk 9 jk ( x )£i& < 1}, where g jk is the 
cometric, i.e., (g^ix)) -1 = (gjk(x)). Note that (|2.5|) with A being the identity operator 
is the sharp Weyl formula (pQ, [12], [9]), and the proof of the more general case just 
follows from a straightforward modifications of that of this special case. 

As a first step in the proof of Theorem l2.1[ let us use some ideas from the proof of the 
Duistermaat-Guillemin theorem [6] (see also [8]). Given (x,^) e S*M we define L(x,£) 
for (x, £) € V to be the minimal t > so that $t(a:, £) = (x, £) and we define £) to be 
+00 if (x,£) ^ P, where P is as in (|2.1j) . Then L(x,£) is clearly a lower semicontinuous 
function on S*M . As a result, since we are assuming that V has measure zero, it follows 
that for a given T > 

P r = {(i,f)6S*Af: i(a;,e)<r} 

is a closed subset of S*M which is of measure zero since Vt C V . Therefore, given e > 0, 
we can find a pseudodifferential operator & <E \&°,(M) whose principal symbol satisfies 
< b (x,£) < 1, b (x,£) = 1 for € Af(P T ), where A/"("P T ) is a neighborhood of 77 

in S*M and 

/ |& (a;,e)| 2 dT/^<e/3. 

JB'M 

By Lemma \2. 21 we conclude from this that 

(2-6) ]T I \bexfdV Ks/S + O.iX- 1 ), 

since we are assuming that |T*73| = 1. 

If we let B = Id — b € \&°,(M) then we claim that there is a uniform constant C, which 
is independent of e and T above so that 

(2.7) sup/ iTJeAj | 2 dV < C/T + C' B T A -1 ^ 2 . 

If T is chosen large enough so that C/T < e/3, the preceding inequalities imply that 
there is an A = A (e) so that 

(2.8) -i- £ SU P / |e Aj ( 2; )| 2 ^<e, ifA>A . 
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As we shall see, this and (|2.4[) immediately yield Theorem 12. II 

Our main estimate (|2.7I) would follow from showing that there is a constant C as above 
so that 

(2.9) sup f \Be x \ 2 ds < CT^A 1 / 2 + C' TB . 

yen J 1 

Here ds is the geodesic arclength measure on 7. Estimate (|2.9j) yields (|2.10[) due to the 
simple fact that if / > then for any 70 G 77 

/ fdV < CX- 1 ' 2 sup [ fds 

for a uniform constant C since 7^-1/2(70) is a tube of width A -1 / 2 about 70. Let 77 c ; 
be the set of unit geodesies that are part of a periodic geodesic. In [20] one of us showed 
that 

[ \e x \ 2 ds = o(X 1/2 ), if 7 eiT\/7 c; , 

which was an o-improvement of the restriction bounds in 4j. The proof of (|2.9[) is an 
adaptation of the one used to establish this result. 

To prove (|2.9p , let us fix a real- valued even function y G S(M) with y(0) = 1 and 
x(t) = 0, \t\ > 1/4, where x denotes the Fourier transform of y. We then have that 

X (T(P-A))e A = e A 

if P = A g . Therefore, in order to prove (|2.9p . it suffices to show that 

(2.10) / l^xm^-A^/l^s^C-T-^^ll/ll^^+C^II/ll^^, j en, 

J J 

where C (but not C' T B ) is a uniform constant independent of T and -B. We shall assume 
in what follows that T is fixed but large, in particular T > 10. 

Note that U is compact. Therefore, in order to prove (|2.10j) . it suffices to show that 
given 70 G JJ there is a neighborhood A/"(7o) of 70 in 77 on which the analog of (|2.10p 
holds with constants independent of 7 £ A/"(7o)- Different arguments are needed for the 
cases where 70 is or is not part of a periodic geodesic of period < T, where T is as above. 

Given 7 G 77 we let T*7 C T*M and 5*7 C S*M be the cotangent and unit cotangent 
bundles over 7, respectively. Thus, if (x, £) G T*7 then £j is a tangent vector to 7 at a; if 
T*M 9 £ — > £fj G TM is the standard musical isomorphism, which, in local coordinates, 
sends £ = G T* X M to £ B - (£j,£jf) with ^ = Efc^C^fc- Note that if 7 G 77 ci 

then L(x,£) = t(j) < 00 for (ac,$) G S**7. With this in mind, we shall let 77 ci (T) denote 
those 7 G 77 cZ for which L(x, < T if (x, ^) G 5*7. 



^Note that, in R 2 , the integral of / > over an 1 X A -1 / 2 rectangle is dominated by A -1 / 2 times the 
supremum of integrals over the line segments in the rectangle that are parallel to the center segment, 
and a similar argument works for the above tubes if one uses Fermi normal coordinates about a geodesic 
which intersects 70 orthogonally. 
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Let us first see that a stronger version of (|2.9p must be valid whenever 7 G 27 C ;(T). 
We first note that if A € *° ; (M) then 

(2.11) A X (T(P - X))f(x) = T- 1 J x(t/T) e~ M ( Ae ltp /) (a) dt, 

and recall that because of the support properties of x(i), the integral vanishes when 
\t\ > T/2. The operator 

f ^ (Ae itp f){x) 
is a Fourier integral operator with wave front set 

(2.12) { (x, t, C, r, y, -77) : $ t (x, = (y, rj), ±r = 0, (a:, G supp A(x, } , 



where p{x,£) is the principal symbol of P = ^—A g and A(x, £) is the symbol of A. If 
i? 7 denotes the restriction to 7 then we are really concerned with the operator 

(2.13) / -> i? 7 Ae l * p /. 

Regarded as an operator from C*°°(Af) -> C°°(7o x [-T/2, T/2]), if supp A(x, f) n 5*7 = 
0, this is a Fourier integral operator of order zero which is locally a canonical graph q If 
7 = 70 G n c i(T) and we take A = B, where B is as above, then this is automatically the 
case since B = Id — b and b has a symbol which equals one in a neighborhood of 5*70 
if 70 G n c i(T). Therefore, by Hormander's [TU] T 2 -estimates for nondegenerate Fourier 
integral operators we have 

T/4 , 

/ \Be* tp f\ 2 dsdt<C\\f\\ 2 L2{M) . 

-T/4 J jo 

The constant C here of course depends on T and 70 (with its main dependence being 
on dist(S'*7o, supp B(x,£)). Since the Fourier integral (|2.13[) with A — B will also be 
nondegenerate if 7 is close to 70, we conclude that whenever 70 G n c i(T), there must be 
a neighborhood 7V(7o) in Pl and a constant C 1o ^b,t so that 

[ T/4 [ \Be up f\ 2 d S dt<C, ^ T \\f\\h { M), 7GAA( 7o ). 

J -T/4 J 7 

If we use the Schwarz inequality and (|2.11|) we conclude from this that 
(2-14) /|i? x (T(P-A))/| 2 d S <TC 70 , B ,T||/l!i 2( Af), 7eA% ), 

which is stronger than (|2.10j) for these 7. 

Let us now see that we also have favorable bounds on 77\77 C ;(T). If we fix a 70 in this 
set and choose a C G ^ Q cl (M) whose symbol vanishes on a conic neighborhood of T*7o 
then by the above arguments there must be a conic neighborhood of 70 on which we have 
the analog of (|2.14l) when B is replaced by Co B. This fact is independent of whether or 
not 70 is periodic. It is just our earlier observation that (|2.13l) is a nondegenerate Fourier 
integral operator when the symbol of A vanishes in a conic neighborhood of T*jq. 



2 Since, for fixed t, e ltp : C°°(M) — > C°°(M) is a nondegenerate Fourier integral operator, one needs 
only to check this assertion for t = 0, in which case it is an easy calculation using any parametrix for the 
half- wave operator. 
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Thus, in order to show that we have uniform bounds as in (|2.10[) on a neighborhood of 
such a 70 G II\II c i(T). it is enough to show that if A G ^° cl (M) has a symbol supported 
in a small neighborhood of T*7o then we have 

(2-15) / \A X (T(P - A))/| 2 ds < CT-'X^WfWl^ + C' T , Aria ||/||| 2(M) 

for every 7 G 77. 

Note that for every x G 70, T x ^$ is one-dimensional and if £ G T*7 then — £ G T*"/q, 
since ±£j G T x j are the corresponding tangent vectors to 70 at x pointing in opposite 
directions. Thus, T*7o naturally splits into two components, which we shall denote by 
TJ.70, and in order to prove (|2.15|) . it suffices to show that the estimate holds if the 
symbol of A is supported in a small neighborhood of one of them, say, Tf 70, since the 
same argument will apply to T*7 . 

We shall assume in what follows that the injectivity radius of (M, g) is 10 or more. If 
not than we can subdivide 7 into a finite number of segments of length smaller than one 
tenth of the injectivity radius and use the argument that follows to prove the analog of 
(|2.15l) for each of these, which in turn yields (I2.15[) for all of 7. 

Let Sf = Ax{T(P - A))/| 7 then we wish to show that 

(II<S1Il 2 (m)^l 2 (7)) < CT^A 1 / 2 + G'T,A,~f - 

This is equivalent to saying that the dual operator S* : L 2 {^) — > L 2 (M) with the same 
norm, and since 

\\S*9\\h(M) = / SS *9Vds < \\SS*g\\ L 2 {l) \\g\\ L 2 h) , 
we would be done if we could show that 

(2-16) \\SS*g\\ L 2 (l) < (cT- 1 A 1 / 2 + Cr,A, 70 ) IMU'W 

But the kernel of SS* is K("f(s), 7(s')) where "f(s) parameterizes 7 by arclength and 
K(x, y), x,y G M is the kernel of the operator A o p(T(P - A)) o A* with p = (x(i~)) 2 
being the square of %. Its Fourier transform p is the convolution of x with itself and thus 
p{t) = 0, \t\ > 1/2. Consequently, we can write 

(2.17) A o p{T(P - A)) oA*=T- 1 [ / p{t/T) e~ ltx (a o e ltp o A*) dt. 
The wave front set of the kernel of 

A o e ltp o A* 

regarded as an operator from C°° (M) to C°° (M x R) is contained in 

(2.18) {(x, t, £, r; y, -77) : $ t (x, = (y, 77), t = ±p(y, 77), (x, £), (y, rj) G supp A}. 

Our assumption that 70 ^ Pl c i{T) implies that if (£,£) G 5'*7o then {$t(x, ^) : 1 < 
\t\ < T} must be a closed subset of S*M which is disjoint from {(x, £)}. If we assume 
also that (x, £,) and (y, 77) belong to the same component 5*170 of 5*70 then we have that 
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®t ( x '0 — {ViV) for some \t \ < 1, and therefore $ t (a;,^) ^ (y,i]), for \t\ E [2,T— 1] since 
if = (y,rj) then we must also have &t-t (x, — ( x i£) f° r this to. Consequently, 

{(x, $ t (x, o) : (x, g s; 7o , 2 < |t| < t — 1} n s; 7o x s; 7o - 0, 

and since both are compact subsets of S*M x S*M, we deduce from (|2. 18|1 that if the 
symbol of A is supported in a small conic neighborhood of 1S+70, then K(t, x, y) will be 
C°° when |t| G [2.T-1]. 

Therefore, for such A, if if a S Cq°(R) equals one if |t| < 3 and zero for \t\ > 4, the 
difference between the kernel K(x,y) in (|2.17[) and 

tf (z, y) = J 1 " 1 /" a {t)p{t/T)e- ltx (A o o A*) (z, y) dt, 

must be bounded, by a constant which is independent of x and y (but depends on T, 
70 and A). Since we are assuming that the injectivity radius of (M,g) is 10 or more 
one can use the Hadamard parametrix construction for the wave equation and standard 
stationary phase arguments (cf. Chapter 5 in [TO] or the proof of Lemma 4.1 in [3]) to 
see that for x, y <E M we have 

\K Q (x,y)\ < CT^X 1 / 2 (di S t g (x,y) )-^ 2 + C A . 

Since this kernel restricted to 7 x 7 gives rise to an integral operator satisfying the 
estimates in (|2.16p . we conclude that we also have uniform bounds of the form (|2.10[) . 
when A is as above. 

This completes the proof that the analog of (|2.10l) holds for all 7 in some neighborhood 
of 70 when 7o G II\II cl (T). 

Combining what we have done for II c i(T) and 77\77 c i(T), since 77 is compact, we con- 
clude that (|2.10[) must be valid with uniform constants for every 7 G 77. This completes 
the proof of (|2.9p and hence (12.8[) . Since the latter holds for all e > 0, we conclude from 
ED that 



(2.i9) \ m ^ p ivTAj £ (\ Tl/8 IK-IUW 12 = o- 



A,<A 



We can now finish the proof of Theorem 12.11 using a counting argument from [5] and 
6]. If S C {Aj}, we define its density to be 



D(S) = liminf 



#{Aj G S : Aj ■ < A} 



a^oc AT(A) 

Then if we use (12 . 19[) we conclude that for every n = 2, 3, . . . we can find a subset S n of 
the eigenvalues so that 

D(S n ) > 1 and A7 1/8 ||e A , ||i*(M) < Aj G 5 n . 

n J n 

Using this we conclude that there must be a set Soo — {Xj k } C {A.,} of density 1 so that 

limsupA ifc 1/8 ||e Aj .J| L 4 (M) =0. 

Indeed, by the above, we can choose increasing N v G N, v = 2, 3, . . . so that 
#{Aj G S v : Xj < X}/N(X) > 1 - 2/v, VA > 7V n _i. 
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Consequently, 

oo 

Soo = IJ s v n {\j < n v } 

will have the desired properties. □ 

3. Average L 4 norms of spherical harmonics. 

In this section we shall prove (|1.8p : 

Theorem 3.1. Let {Y^(x)}t,__ k , x = (sin cos #, sin </> sin 6*, cos <?!>), be the orthonormal 
basis of spherical harmonics of degree k defined in (jTTTJ) . Then there is a uniform constant 
C so that 

(3-1) 2fcTT ^ 1 2 \Y„ k fdV <C\ogk, k>2. 

m— — k 

Moreover, if 11 = (0, 0, 1) and if r = min± dist(x, ±1) 

k 1/4 fC^/V-^iog^r)) 174 , r>2/k 



(3-2) ( J2 l^)| 4 ) 1/4 < 



Ck 1 ' 2 , r<2/k. 



Clearly (|3.2p implies (|3.ip , and so we just need to prove the second inequality in the 
theorem. To prove this, we first realize that by Parseval's theorem we have 

d6». 



2^ £ |^(x)| 4 = /I £ 



The kernel LJ k (x, y) for projection on to spherical harmonics of degree k is given my the 
formula 

fe 

n k (x,y)= Y*(x)Yj(yj, 

m— — k 

which means that 

(3-3) 2tt J2 \ Y ™(*)\ 4 = \n k {x,e^x)\ 2 d6, 

rn — — k 

if we abuse notation a bit and let e x denote rotation of our vector x = (xi,X2,xs) by 
angle 9 about the x 3 -axis, i.e., e l9 x = (cos 6x\, sin0x 2 , x 3 ). Using the well known bounds 
(see [33], QI]) for n k , 

\n k (x,y)\ < Ck^ik- 1 + di S t(x,y))-^ 2 , x,yeS 2 

we conclude that 

(3.4) l^nWI 4 <Ck (k- 1 +dist(x,e l9 x))- 1 d8. 

Since dist(a;,e ) < Cr\ sin 6*| , we conclude that the right side of (|3.4[) is < Ckr^ 1 log(fcr) 
if r > 2/k, and < Ck if r < 2/k, for some uniform constant C when k > 2, which is just 
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We believe this estimate to be sharp, but defer the analysis to the future. 

4. L 4 norms of orthonormal bases of spherical harmonics. The space of ONBk 
of Hcrmitian orthonormal bases oi'Hk may be identified with the unitary group [7(2fc+l). 
Any orthonormal basis $k = {4>i, . ■ . , 02fe+i} can be obtained by applying an element 
U E U(2k + 1) to the standard orthonormal basis {Y,^}. We then consider the functional 
on ONB k defined by, 

2fc+l 

A£(**)= H^mllt- 
m— 1 

We have just proved that Af of the standard orthonormal basis is bounded by k log k. 

One may consider similar functionals on orthonormal bases of all eigenspaces with 
Afc < A, i.e. the direct sum ©j.<^^fe- The functional then has a natural generalization 
to any (M, g) and is essentially the one studied in previous sections. 

5. Random orthonormal bases of spherical harmonics. We now consider the At 
functional on a random basis of spherical harmonics. The question we pose is, what is 
the average value of the functional on random orthonormal bases? In [SZj we considered 
problems of this kind for L p norms of individual eigenfunctions, but there is a new 
dimension to the problem for random orthonormal bases. For background on random 
orthonormal bases we refer to |SZJ . 

We introduce the probability space (ONB, dv), where ONB is the infinite product of 
the sets, and v — Hfe!Li v k, where v\ is Haar probability measure on ONBk- A point of 
ONB is thus a sequence $ = {(4>i, ■ • ■ , </ , 2fc+i)}fe>i °f orthonormal bases of Hk- 

The functionals we are interested in are 

2k+l 

(5.5) A£(*)= £ / \\^(x))\\ 4 dV . 

If we fix the standard ONB Y k = {Y^} and express every other as UkY k , then our 
functional is 

2k+l 

(5.6) At(U) = 

Let d/j,k be normalized Haar measure on U(k) and let denote expectation with 
respect to this measure. We conjecture that 

(5.7) E fe A 4 =2(2fc + l), 
i.e. the elements on average have L A norm equal to 2. 

We briefly sketch the proof. We start from the fact that 

(5.8) 

(Su&k+l) U L U m 2 U L U mJ^k{U)) . 



/ WiUY^fdV. 
Js 2 
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In fact, the sum over j is constant, so the right side equals (2k + 1) times 
E (/ Y^(x)Y k m2 {x)Y^{x)Y k mi (x)dV 



The integrals (ju^k+^U^U^U^Jf^dfikiuf) were first studied by Weingarten 

[W| . The main result is that the random variables {\/2k + lUij} behave asymptotically 
like independent complex Gaussian random variables of mean zero and variance one. 
Exact formulae are given in |CSj , and the latter can be used to determine the asymptotics 
of our sums over 2fc + 1 indices with different coefficients as k — > oo. The dominant terms 
come from the cases where all rrij are equal (then one has the fourth moment of the 
Gaussian) or when the indices m,j are paired into couples (one barred and one unbarred). 
Then we have, 
(5.9) 

(2fc + l) 2 EA^(E/) 

~ (2fc + 1) (2 \Yl{xWdV + 2YL^m 2 =- k Us* \Y£M?W*V)) 

= 2(2fc+l) J s2 \n k (x,x)\ 2 dV = 2{2k + I) 3 . 

Here, we use that the 4th moment of the complex normal Gaussian equals 2 and that 
there are two ways to pair the indices in the off diagonal terms. Dividing by (2k + l) 2 
then implies the result. 

6. Other orthonormal bases. Theorem 13.11 shows that the A| functional on the stan- 
dard basis {Y^} is only logfc higher than for a random orthonormal basis. Hence it is 
doubtful that it does not maximize . We do not know which ONB maximizes the func- 
tional, but in this section we suggest a possible construction of one which has a higher 
A^ value than the standard basis. 

As mentioned above, the highest wight spherical harmonic Y£ has L\ equal to k 1 ! 2 
on S 2 , thus maximizing the norm functional. This suggests constructing orthonormal 
bases consisting in part of highest weight spherical harmonics with respect to a well- 
separated set of closed geodesies 7^. That is, for each closed geodesic 7, one introduces 
the subgroup G 7 of rotations fixing 7 (as a set) and then constructs Y^'s with respect 
to this circle action. 

Of course, the are not orthogonal, and their inner products (0*.,$^) depend on 
the angle dj.k between the geodesies 7j,7i- To construct an orthonormal basis it would 
be necessary to apply Gram-Schmidt to such <j)^. , and in the process one may destroy 
the high L A norms of the resulting eigenfunctions. The question is, how many <j>* can be 
used in such a construction while preserving the high L 4 norms of these Gaussian beams? 

The geodesies 7^ are points in the space G(S 2 ,go) of geodesies of S 2 . A well-separated 
set of 2k + 1 geodesies (i.e. a basis) would only have separation of order To 
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beat the bound k log k for the standard basis one it would suffice to construct a partial 
orthonormal basis containing k x ~ 5 roughly Gaussian beams with roughly ||(/>^ \\\ ~ \fk 
and with S < \ . One would then complete it with an arbitrary orthonormal basis of the 
ortho-complement of the span. It would be interesting to see how far separated the jj 
would need to be so that Gram-Schmidt would not destroy the bounds \ \\ ~ \fk too 
much. 
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